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Abstract 

We prove that the scalar and 2x2 matrix differential operators 
which preserve the simplest scalar and vector-valued polynomial mod- 
ules in two variables have a fundamental Lie algebraic structure. Our 
approach is based on a general graphical method which does not re- 
quire the modules to be irreducible under the action of the correspond- 
ing Lie (super)algebra. This method can be generalized to modules 
of polynomials in an arbitrary number of variables. We give generic 
examples of partially solvable differential operators which are not Lie 
algebraic. We show that certain vector-valued modules give rise to 
new realizations of finite-dimensional Lie superalgebras by first-order 
differential operators. 
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1 Introduction 



There has been a significant amount of interest recently in the study of 
Schrodinger operators which preserve exphcit finite-dimensional subspaces 
N of the underlying space of smooth functions, p5|, 0, |[Tl|]. A Schrodinger 



operator, or more generally any differential operator T, which enjoys this 
property is said to be partially solvable. Indeed, if M is n-dimensional and 
the restriction of T to A/" is a self-adjoint operator with respect to a suitable 
inner product, one can then determine n eigenvalues and eigenfunctions of T 
algebraically, counting multiplicities. Notable examples of partially solvable 
Schrodinger operators in one dimension include some new families of anhar- 
monic oscillator potentials, as well as the familiar Morse, Poschl- Teller and 
Mathieu potentials, which have been thoroughly analyzed from this algebraic 
point of view, [^, |jT9|, [|12|, 0]. It has been observed that most of the known 



examples of partially solvable Schrodinger operators can be written as ele- 
ments of the universal enveloping algebra of a finite-dimensional Lie algebra 
of first-order differential operators admitting a finite-dimensional subspace 
of the underlying Hilbert space as a module of smooth functions. The dif- 
ferential operators satisfying this property are called quasi- exactly solvable, 
or more generally Lie algebraic if no assumption is made on the existence of 
such a finite-dimensional module. One can thus construct many examples of 
partially solvable Schrodinger operators by first classifying the Lie algebras 
of first-order differential operators admitting finite-dimensional modules of 
functions, and then looking for elements in the universal enveloping algebra 
which can be put in the Schrodinger form after conjugation by a suitable 
non- vanishing multiplication operator and a local change of the independent 
variable. This conjugation is referred to usually as a gauge transformation. 
Note that the gauge transformation will generally not be unitary, so that 
square integrability need not hold at the outset, but only for the eigenfunc- 
tions of the resulting Schrodinger operator obtained after conjugation. This 
program has proved to be remarkably successful in that it has enabled one to 
construct many new examples of quasi-exactly solvable potentials in higher- 



dimensions, pO[, 0, as well as matrix potentials for particles with spin, |jT9 
i, i- 

For all the examples known so far, the invariant modules have a partic- 
ularly simple structure when expressed in a suitable coordinate system and 
gauge. In the one- dimensional case, they are the vector spaces of polynomi- 
als of degree less or equal than a positive integer n, and the corresponding 



2 



Lie algebra is simply the standard representation of sl2 by first-order differ- 
ential operators in one variable. It is easy to show that there are no other 



possibilities in one dimension, |T3|, |T^. In the case of quasi-exactly solvable 



differential operators in two variables associated to a complex Lie algebra, 
the modules also take the form of polynomial modules in the appropriate 
coordinates and gauge, 0. They are either one of the vector spaces of bi- 
variate polynomials x^y^ of bidegree less or equal than an ordered pair 
(n, m) of positive integers, or one of the vector spaces of polynomials of total 
degree i+j less or equal than a positive integer n, or one of the vector spaces 
of polynomials whose bidegree satisfies a linear constraint of the form 
i + rj < p, for some positive integers r and p. We will refer to these modules 
as the rectangular, triangular or staircase modules in view of their structure 
within the lattice of bidegrees. The corresponding Lie algebras are given by 
the standard realizations of ©£[2, sfs and ^[2 k C"*"^ by first-order differen- 
tial operators in the plane. The inequivalent real forms of these algebras give 
rise to different quasi-exactly solvable potentials, [|10[. It should be noted 



that while the rectangular and triangular modules are acted on irreducibly 
by 5I2 © ^[2 and sl^, this is not the case for the staircase module. 

A fundamental question is to determine why practically every known ex- 
ample of partially solvable differential operator should be expressible as a 
polynomial in the generators of a finite-dimensional Lie algebra of first-order 
differential operators leaving the corresponding subspace invariant, and to 
what extent such a polynomial is unique. In particular, one would like to 
know if there exist modules for which the generic partially solvable differ- 
ential operators are not Lie algebraic. These are the main problems that 
we address in this paper for the case of rectangular, triangular and staircase 
scalar or rank 2 modules. In the case in which the module is acted on irre- 
ducibly by the Lie algebra, which in our case would be the rectangular and 
triangular modules, a partial answer is provided by the Burnside Theorem, 



as pointed out by Turbiner, . Recall that according to the Burnside The- 
orem, if y is a complex vector space and a is a subalgebra of End V acting 
irreducibly, then every endomorphism of V can be expressed as a polynomial 
in the generators of a. This implies that any differential operator preserv- 
ing a rectangular or a triangular module can be expressed as the sum of a 
polynomial in the generators of the realization of the corresponding Lie al- 
gebra and a differential operator annihilating the module. Note, however, 
that this argument does not provide any direct information on the relation 
between the order of a differential operator leaving a rectangular or triangu- 
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lar module invariant and the minimum degree of a polynomial representation 
thereof, owing to the presence of polynomial relations between the generators 
of the Lie algebra, which form a basis for the primitive ideal associated to the 
module. The case of the staircase modules is considerably more complicated, 
since the underlying Lie algebras of first-order differential operators do not 
act irreducibly, so that one cannot apply the Burnside Theorem. In fact, 
we will see that in stark contrast with the rectangular and triangular cases, 
the generic differential operator preserving the staircase module is not Lie- 
algebraic. One needs therefore to develop a different approach altogether. 
This is what we suceed to do in this paper, by developing a general graphical 
method which does not require the (scalar or rank 2) modules to be irre- 
ducible under the action of the corresponding algebra. We will thus entirely 
bypass the use of the Burnside Theorem and the need for extensive calcula- 
tions in local coordinates. Indeed, the latter can get prohibitively difficult in 
higher dimensions. (Even the one-dimensional scalar case is quite involved 



when one works in coordinates, [^). It is worth noting that the graphical 
method can be generalized without major difficulties to the case of rank A; 
modules of polynomials in variables. It is thus a powerful tool for future 
applications. 

In Section 2, we briefly review the normal forms for the maximal Lie 
algebras of first-order differential operators in two complex variables which 
admit finite-dimensional modules of polynomials when expressed in a suit- 
able coordinate system and gauge. In Section 3, we unravel by means of a 
general graphical method the Lie algebraic structure of the scalar differential 
operators of arbitrary finite order which preserve a rectangular, triangular 
or staircase module. We give a formula for the minimum degree of a polyno- 
mial in the generators of the corresponding Lie algebra representing the given 
quasi-exactly solvable differential operator. In the case of triangular modules, 
we also obtain an exact formula for the number of free parameters which de- 
termine the most general partially solvable differential operator of arbitrary 
finite order. Section 4 is concerned with 2x2 matrix differential operators of 
arbitrary finite order admitting the direct sum of two rectangular, triangular 
or staircase modules as an invariant subspace. For suitable values of the 
gap parameters, the underlying algebras of matrix differential operators give 
rise to new realizations of finite-dimensional Lie superalgebras of first-order 
differential operators. The abstract structure of these finite-dimensional Lie 
superalgebras is identified in every case. For the remaining values of the gap 
parameters, one obtains an infinite-dimensional Lie superalgebra. Remark- 
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ably, the graphical method introduced for scalar modules enables us in both 
the finite and infinite-dimensional cases to prove structure theorems for the 
corresponding partially solvable matrix differential operators. 

2 Quasi-exactly solvable Lie algebras in 

Our purpose in this section is to briefly recall from ^ the normal forms for 
the maximal Lie algebras of flrst-order differential operators in two complex 
variables admitting flnite-dimensional modules of polynomials in suitable co- 
ordinates, along with the structure of these modules. 

Let U be an open subset of C^, with local coordinates {x,y). The flrst- 
order differential operators on U, 

T = fix, y)d^ + g{x, y)dy + h{x, y) , 

with coefficients in C°°(f/, C) form a Lie algebra which we will denote by 
25(1) _ 2)(i)(f/). The space C°°([/,C) is naturally a S)(^)-module. We will 
use the terminology of 0, and refer to finite-dimensional subalgebras g of 
D^^^ admitting a finite-dimensional module J\f C C°°{U,C) as quasi-exactly 
solvable Lie algebras. The local diffeomorphisms ip of U and the rescalings 
by smooth functions e°" G C°°(f/, C*) define an infinite pseudogroup G of 
automorphisms $ of D^^\ whose action is given by 

$(T) = ■ (^,T ■ e"" . (2.1) 

The classification under the action of G of all quasi-exactly solvable Lie 
algebras in two complex variables is known, |^ (its real counterpart has been 
recently completed, 0]). We shall focus our attention on the three maximal 
families of equivalence classes of quasi-exactly solvable Lie algebras whose 
normal forms preserve a module of polynomials. 

i) qI^^ ~ © SI2, with generators: 

'-^n 1 5 '^n xOx ~ , 

K+ = y%-my, R- = dy , Kl = ydy--, (2.2) 

with n,m E Z"*", and associated finite-dimensional module given by 
7^n,™ = {x'y^ \ 0<i<n, 0<j <m}. 
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ii) — sis, with generators: 

JI = x^d^ + xydy -nx , = xydn, + ffdy - ny , 

J3=|/a,, J' = a,, f' = dy, J^=Xdy, 

Jl = xd,-^, J^^=ydy-^, (2.3) 

with n G Z+ and corresponding finite-dimensional module given by 

X = {xy I < i + J < n} . 

iii) Qp"^'^ ~ X C^^, with generators: 

= d^, Jp = x^d^ + rxydy - px , = xd^ - -, 

J4 = yQ^ ^ j5+i = < i < r , (2.4) 

with 2 < r G Z+, G Z+ and, for g G Z+, associated module given by 

S;,, = Wy' I < « + rj < p , < J < g} , 

As a special case, we shall consider the module S"^ = S^p^p/^]- 

The nomenclature of the algebras is the same as the one used in ^j, and 
the letters 71, T and S used to denote the modules reflect of course the 
rectangular, triangular and staircase structure of the lattices in the positive 
Z+ X Z+ quadrant corresponding to the powers (i, j) that appear in 7ln,m, 
Tn and S^ ^. The integers n, m, p, q which label the reahzations of SI2 ©-SI2, 
sis and q\.2^ C"*"^ in the list above are related to the global models for these 
algebras, where the underlying complex surfaces are respectively given by 
Pi(C) X Pi(C), P2(C) and the r-th Hirzebruch surface S,., and the modules 
are sections of holomorphic line bundles over these surfaces, j^. 

3 Scalar k-th order differential operators in 
with invariant modules of polynomials 

Our aim in this section is to prove structure theorems for the subalgebras 
of the associative algebra 2) = S)(f/) of scalar differential operators in two 
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complex variables 

T= hi^,y)didi, f,,ix,y)EC^iUX), (3.1) 

which map the modules 7ln,m, %i and ^ to themselves. Let us remark that 
our results will also be valid when x and y are restricted to the real plane M^. 
The motivation for our study arises from a theorem of Turbiner concerning 
scalar k-th order differential operators in one variable, To set Turbiner's 
theorem in context, let us recall, [T^, [Q, that up to equivalence under the 
one-dimensional analog to the pseudogroup G defined in ( |2.1] ), there is a 
unique family of quasi-exactly solvable Lie algebras of first-order differential 
operators on the line, namely the family of representations of given by 

fln = {Jn, J~: Jn]^ ^hcrC 

J+ = x'^d^ -nx, J- = d^, Jl = xd^ - -. (3.2) 

The associated invariant module is the polynomial module 

Vn = {x'\Q<i<n}. 

Let D be the associative algebra of all linear differential operators in one 
complex variable, and let S*^*^-* denote its subspace of differential operators of 
order at most k. We shall denote by the subalgebra of J) spanned by all 
differential operators of the form T(9^~^^, where T is an arbitrary differential 
operator. Let be the subalgebra of D of differential operators which 
preserve and let ^pi''^ = ^„ n Finally, let $^{s{2) = ^Lo^Ksk) 

be the direct sum of the first k symmetric powers of s[2. Turbiner's theorem 
may be reformulated as follows: 

Theorem 3.1 i) The linear map : S'^(s[2) '^n^ determined by the 
representation pn : SI2 — > 0n is surjective for all k = 0, . . . ,n. ii) If k > n, 
then = %i^n^ © (!)('=) n SO")). 

Remark. Turbiner's original proof was done by explicit computation. We 
present here a simplified proof whose principle will serve as a guide in the 
higher dimensional case. 

Proof. We first note that kervr^ = J)*^^"-*, where 7r„ denotes the induced 
homomorphism 7r„ : ^„ —>■ End('P„), so ii) clearly holds. Therefore, the 
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restriction of 7r„ to ^i"''' is injective, so dim^n"'' < (^ + 1)^- Now, the mono- 
mials {(^n )*(^°)-'~*}i=0' J = 0) 1) • • • ) all linearly independent differential 
operators. In particular, there are (n + 1)^ monomials of order at most n, so 
they form a basis of ^i"''. Since = x^^'^dl + ■ ■ ■ , no monomial of 

degree higher than k may appear in the expansion of a differential operator 
T^'^) e Q.E.D. 

Thus, Theorem ^TT| simply expresses the fact that any differential operator 

k<n 

T = £/,(x)9:, 

i=0 

preserving Vn may be written as a polynomial of degree k in the genera- 
tors CO). 



Later, Turbiner proposed a proof of Theorem based on the irreducibil- 
ity of the module Vn under the action of g„ and the Burnside Theorem. 
However, this argument only guarantees that T'^'^^ will be the sum of a poly- 
nomial in the generators of Qn with an operator annihilating the module. In 
particular, it gives no information on the degree of the polynomial. 

The linear map : S^{sl2) ^^f"* C ©'■'^^ is clearly not injective for 
k>2. In fact, the set )*(-'^)-'~*}i=0' J = 0; • • • , ^ is a basis of imp^ (the 
completeness follows from the quadratic relation J~ = JnJn ~ Jn~ \n{n + 
2) which expresses the scalar action of the Casimir operator on Vn)- Thus, 

dim(kerp^) = dim§'^(s[2) - {k + if = hk + l)k{k - 1) . (3.3) 

6 

The dimension of the kernel of the induced map p^ : 6*^(312) 
may also be easily obtained: 

dim(kerp^) = ^A;(A;- 1). (3.4) 



We emphasize that equations (3^) and (|3.4| ) are still valid if the cohomology 



parameter n labeling the representation 0„ is replaced by an arbitrary com- 



plex number A; see [15] and M for more details. 



We now focus on the two-variable modules TZn,m, %i and S^ ^, and define 
^n,m, and 6p g as the subalgebras of D which preserve these modules. Let 
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T^n,m ■ 5^n,m ^ End(7e„,m), 7r„ : T„ ^ End(7;), and vr^ ^ : Q^ ,^ End(5^^g) 
denote the induced homomorphisms of associative algebras. We have the 
following elementary result: 

Lemma 3.2 The kernels of 'Kn,m, 7r„, and vr^g are given by 
ker7r„,„ = {T G D | T = ^ f,,{x,y) dld^} , 

i>n or 
J >m 

kervr. = {T E ^ \ T = S^S^} , 

i+j>n 

ker7r;, = {TGS)|T= ^^S^} . 

i+rj>p 
or j>g 

We now choose the following distinguished complements to the kernels in 
Lemma p.2| : 

0<i<n 
0<j<m 

T„ = {TGTJT= Yl Ux,y)dldl}, 

G;^^ = {Te&;jT= J2 U^,y)dldi}. 

0<i+rj<p 
0<j<q 

Note that the order of a differential operator in the distinguished complement 
&pg cannot exceed p. The next result follows from a simple constructive 
argument: 

Lemma 3.3 The restrictions Tin^m '■ 9^n,m End(7?,„^m)) '■ '^n ^ End(7^), 
and vfp g : (5^ ^ — > End(5p ^) are vector space isomorphisms. 

Proof. Let us prove the lemma for one of the modules, say 7^. Let 7 G 
End(7^). The differential operator T G X„, with coefficients fij given by 

/oo = , /lo = T(x) - x/oo , /oi = 7{y) - yfm , . . . 

clearly satisfies TCniT) = 7. Q.E.D. 
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Note that the coefficients fij of the differential operators in the distin- 
guished complements are polynomials in x and y. The subalgebra Dp oi 
D of differential operators ( p. ID with polynomial coefficients fij inherits a 
natural Z x Z grading from C[x,?/]. We shall say that a differential oper- 
ator T G Dp has bidegree degT = if T{x'°y^°) G {x^'^+'y^o+j'^ foj. 
(^O)Jo) G Z+ X Z+. Note that the generators of qI^ g^'*''' all have 
well-defined bidegree, [^. A set of differential operators in Dp with different 



bidegrees is of course linearly independent. 

We will obtain in what follows structure theorems for each of the distin- 
guished complements. Let us remark that both q}^^ and q}^ act irreducibly 
on their associated modules 7ln,m and 7^. Therefore, a direct application 
of the Burnside Theorem shows that any differential operator T in '!Kn,m 
(respectively Tn) may be constructed as a polynomial in the generators of 
qI^^ (respectively q}^) plus an element of ker TTn^m (respectively ker7r„) 



Note that the action of Qp'^''^ on Sp ^ is neither irreducible nor completely 
reducible, contrary to the the assertion in ref. ||2^. In fact, we will exhibit 
differential operators in &p ^ which cannot be expressed as a polynomial in 
the generators of q^''^ plus an element in ker vr^ g. 

3.1 The rectangular module IZn^m 

Let denote the subspace of D spanned by all differential operators of 

x-order at most k and y-oidei at most /, and let '^n-m = ^n,m H 2) 
The structure theorem for !H„,m is a straightforward generalization of Theo- 
rem 



3.1 



Theorem 3.4 Let pHx) : ^^(sfc) ^„(x) and pi{y) : S^ish) ^ ^™(l/) he 
the linear maps determined by pn{x) : SI2 — > Qn{x) and Pm{y) '■ sfc — > Qm{y), 
respectively. Then, 

0<i<A; 0<i<fc 

o<i</ o<j<i 
is surjective for all k = 0, . . . ,n and 1 = 0,... ,m. 

3.2 The triangular module 7^ 

Let D^'^^ be the subspace of D spanned by all differential operators of order at 
most fc, and let T^f'^ = f„nD('=). As before, let S''{sh) = 0j=^o S'{sk) be the 
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direct sum of the first k symmetric powers oish, and let : S^s^) 
be the hnear map determined by the representation p„ : sis dl^- Our 
purpose is to construct a suitable basis of impj^ for all k G Z+, from which 
the structure theorem for T„ will arise as a simple corollary. The bidegree of 
the generators (|2.3|) of q!^ is given by: 

deg = (1, 0) , deg 4 = (0, 1) , deg = (-1, 1) , deg = (-1, 0) , 
degJ^ = (0,-1), deg J6 = (1,-1), degj: = (0,0), deg J„« = (0, 0) , 

and the bidegree of a monomial = {JnY^ ■ ■ ■ (Jn)''^ is simply Y^^=o deg Jn 
Let \K\ = ki + ■ ■ ■ + kg denote the degree of . Note that \K\ is just the 
order of . We define the length of a monomial as 

I J^l = min \L\ , 

deg J^=deg 

where the minimum is taken over the set of monomials of degree at most \K\. 
In other words, if deg = then | is just the minimum number of 

"steps" required to map a point (zq, jo) ^ Z+ x Z+ to {io + i, jo + j)- A mono- 
mial is of maximal length if | = l-ft'l- The monomials 1, J^, . . . , are 
obviously of maximal length. The following facts follow from an easy graph- 
ical argument, cf. Fig. 1: 

i) There is exactly one monomial of maximal length of bidegree for 
each G Z x Z, up to the ordering of its factors. 

ii) The monomials of maximal length / are {{Jn^^)-' (^^)'~"'}j=i_5,(,, with 
s = 1,... ,6 and jo = J6. 

• • • 

• o o • ® Length / = 

• o ® o • o Length I = 1 

(0,0) 

• o o • • Length I = 2 



Figure 1: Lattice of bidegrees of monomials of length / < 2. 

Note that there are 6/ monomials of maximal length I > 1. If is a 
monomial of maximal length, the monomials Jn {JjiY {JnY ^ ^'j — 0' 
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linearly independent, and have the same bidegree as J^. Remarkably, the 
monomials of this form are also a generating set of imp^ for all k G Z"*", 
according to the following lemma: 

Lemma 3.5 

i) The monomials 

{iJn-yiJny-'i4)VnY''-' \ 1 - ^10 < J < I < t < k , < Z < t ^ l} , 

(3.5) 

with s = 1, . . . ,6 and J° = J^, form a basis of im for all k. 

a) The number of k-th order monomials in the basis is Uk = {k + 1)^, and 

dim(imp^) = ^(fc + l)2(fc + 2)2. (3.6) 

Proof. Linear independence of the monomials follows from the preceeding re- 
marks. Completeness is a consequence of the following 9 quadratic relations, 

/ tT \ 2 tT t8 t7 ^ t8 ^ / ^ \ 
\Jn) "I" Jn^ri ~ '^n ~^ S" ~ "S "S ' 

-^(^ + 1), Jlj' = j'n{4 + 1), 

(3.7) 

which allow us to reduce any monomial to a linear combination of the el- 
ements of the set (|3.5| ). The second assertion now follows by straightforward 
computation. Q.E.D. 











= J (j„ + ~ ("^ 


+ 1)), 






f - 










= J'{JI + 1 + 1). 





From Lemma p.5| , we obtain the following result: 

Theorem 3.6 The linear map p^ : S'^(3[3) — ^ determined by the repre- 
sentation pn : sl^ g}^ is surjective for all k = 0, . . . ,n. 

Proof. The map p^ : S"(s[3) T„ is clearly surjective, since dim'T^ = 
+ l)'^{n + 2)^. The highest order terms of the monomials in (|3.5| ) are all 
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independent, which imphes that no monomial of degree higher than k may 
appear in the expansion of a differential operator T^^) e Q.E.D. 



Theorem 3.6 is stated in a different form in [23|. It should be noted, 



however, that just as in the SI2 case, this theorem does not follow directly 



from the Burnside Theorem, contrary to what is stated in |23|. As another 
immediate consequence of Lemma |3.5| , we have: 

Corollary 3.7 The kernel of the homomorphism pn : ^(sls) T) deter- 
mined by the representation pn : sta — ^ 0^^ is the primitive ideal generated by 
the quadratic relations p.7| ). 

Let us emphasize that these results are completely independent of the 
discrete character of the cohomology parameter n labeling the Lie algebras 
qI^, being also valid for a complex cohomology parameter A, as considered 
in 0] and The graphical method used to obtain the basis ( p.5|) can be 



generalized to the A^-dimensional representations of sIn+i which appear as 
the hidden symmetry algebra of the Calogero model and its quasi-exactly 
solvable extensions, 0, [|T^, Finally, let us give the formula analogous 
to (pD for the map : S^isls) dW/D^^-i): 

dim(kerp^) = i^l±lM^lZ_ll (/^4 ^ + 323k^ + 1988A; + 2052) . 
3.3 The staircase module SI,„ 

The structure of the space of operators preserving this type of modules is 
considerably more complicated than that of the other two modules. As re- 
marked before, the action of q^'"^ on S^^^ is reducible, leading to a number of 
interesting consequences. 

The bidegrees of the generators of g^^'^ are: 

deg = (-1, 0) , deg Jl = (1, 0) , deg Jl = (0, 0) , 
degJ^=(0,0), deg J^+* = (i, -1) , 0<i<r. 

We shall refer to the second component of the bidegree as the ^/-degree. Note 
that the y-degree of the generators (|2.4| ) of q^'"^ is nonpositive. Thus, no 



operator in the distinguished complement (3^^ with positive y-degree (as 
T = yd^, if r = 2, p > 2 and q = [p/2]), can be obtained as a polynomial 
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in the generators plus an operator in kervr^^. Moreover, not every operator 

in &pq with nonpositive y-degree may be written as a polynomial in the 
generators of 0^^'^" only. An example is given by T = a;^(x^(9^— 2x9^— 2^(9^+2) 
in the case p = 2q = r = 2. Our next objective is to prove that every k-th 
order differential operator T^^^ in ©p^g, the subspace of ^ of differential 
operators with nonpositive y-degree, may be written as the sum of a k-th 
degree polynomial in the generators of flp^''" and an operator in kervTp^. We 
shall construct a suitable collection of monomials in the generators ( p.4j ) 
such that their projections in ^ along ker 11^^^ form a basis of ©p|g- (Note 
that in the triangular case every monomial of degree at most n has a non- 
trivial projection in T„ along ker vr^). We follow the approach of the previous 
section, obtaining first a set of monomials of maximal length (in this case the 
bidegree does not determine uniquely the monomials of maximal length), and 
then extend it to a basis of the image of the map p'^'' : S'^(0t2 ^ C^'"*"^) — > SD^'^^ 
determined by the representation Pp : ix C"*"^ flp^'*^- 

Lemma 3.8 

i) The monomial^ 

{(j;+^)'-^(j'+^7(jp')"(j^)*-'-", (j^)v^+'(j^+'^)'"^-^(jj)"(j^)*-'-"}, 

with e = 0,1, < j < I < t < k, < s < I - 1, < n < t - I, and 
1 < i < r — 6so form a basis of im p^''^ for all k. 

a) The number of k-th order monomials in the basis is 

nk = l{k + l){k + 2) ((r + 2)k + 3) . 


The basis described in the previous lemma does not completely fulfill our 
requirements, for it contains monomials of order less than p acting trivially 
in Spg. A subset of monomials acting non-trivially can be constructed by a 
simple graphical argument, cf. Fig. 2. 

Lemma 3.9 The projections in along ker tt^^ of the monomials 

{{jp^yiJ^-^^yiJpTiJ^, {j'yj'+\j'+y-'-\j^)''{j^y"'}, (s.s) 

^By convention, the second group of monomials is not present if Z = 0. 
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® Length / = 
(0, 0) o Length / = 1 

• o o o • ® • Length I = 2 

® Length I = 3 
• Length / = 4 



Figure 2: Lattice of bidegrees of monomials acting non-trivially in 

with e = 0,1, < j < q, < s < p - jr, < t < J - 1, 1 < i < r - 6to, 
< n, < m < g — j, and n + rm < p — s — jr (where s = s for the first 
group of monomials, and s = for the second group), form a basis of &p^q- 

Proof. The projections of the monomials in 6^ ^ along ker tt^ ^ are certainly 
linearly independent. The number of monomials in (p.8|) is given by 



<? P-jr _ ^ 

= 5Z ^)^{^^P- ^ - "^in [ ]) - j) 

j=0 s=0 



where /(j, s) = 2 if s > 1 and /(j, 0) = jr + 1, and <j(ro,Po, Qo) = dim^^o 
(go + l)(po + 1 — ^f^oQo)- We leave as an exercise to the reader to verify that 



N;^ = dim©;^^ = ii±lKi±^(^(p+i)(p+i_g^) + ^(3g + l)) . Q,E.D. 

The structure theorem for (5^^^ now follows from the fact that the highest 
order derivatives of the monomials (|3.8| ) are all independent. 



Theorem 3.10 Let T^^^ he a k-th order differential operator in &^p\- If 
k < q, then T^^^ may be represented as a k-th degree polynomial in the 
generators of Q^'^ . If k > q, then T^''^ may be expressed as the projection of 
such a polynomial in ^ along kerTr^ ^j. 

We shall now restrict ourselves to the "untruncated" modules Sp = j^^^j 

and the distinguished spaces of differential operators 6^ = (5^ jp^^j which pre- 
serve them. We are interested in obtaining a method for constructing the 
most general differential operator in (5^ irrespective of their y-th degree. Fol- 



lowing |g3[, we introduce a well-adapted grading for the differential operators 
in ©p. If T G 6p has a well-defined bidegree («, j), we define its total degree as 
DegT = i + rj. The constructive method outlined in the proof of Lemma |3^ 
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allows us to obtain the differential operator T representing any given endo- 
morpliism of S^. However, it does not give any information on the order of 
T, which limits its value in quantum mechanical applications. The following 
elementary lemma partially describes the subspace &p^''^ = &p fl D^^^ of k-th 
order differential operators. 

Lemma 3.11 Let T^^^ be an element of Q'p^^^ homogeneous with respect to 
the total degree. Then —kr < DegT*^^) < k. 

Proof. The lower bound is obvious, for T^'^^ is a linear differential operator 
with polynomial coefficients. Now, if DegT*^'^-' = d > k, T^'^^ must annihilate 
the last d diagonals {x'^y^ \ i + rj = p — d + 1, . . . ,p}, but a (nonzero) fc-th 
order differential operator cannot annihilate more than k diagonals. Q.E.D. 

We can use this result and Theorem |3.1CI| to construct the most general 
differential operator in Gp^''^ As an application, we give the explicit form 
of the most general differential operator T^"^^ in &p'^'^\ which may be used to 
construct partially solvable Schrodinger operators in two variables by solving 
the corresponding equivalence problem, . We assume that p > 4 (otherwise 
some of the terms would just have a non-vanishing projection on ker vr^ along 
el). Then, T(2) = J^t.^ T,, where: 

T_4 = aidy, 

T„3 = a2 d^dy + as xd^ , 

T_2 = 04 + 0-5 xd^dy + ae x^dy + a-j ydy + agdy, 

T_i = ag xdl + aio x^d^dy + an yd^dy + ai2 x^d^ + 013 xyd^ + au + ais xdy , 
To = x^dl + ai7 ydl + ai8 x^d^dy + aig xyd^dy + 020 x^dy + 021 x'^ydy 

+ a22 y'^dy + 023 xdx + «24 x'^dy + a25 ydy + 025 , 
Ti = (a27 x'^d^ + 028 ydx + ^29 x^dy + a-^o xydy + a^i x) {xd^ + 2ydy - p) , 
T2 = (032 + 033 y) {xd^ + 2ydy -p+ l){xd^ + 2ydy - p) . 

Note that T_4, ... ,To are linear combinations of monomials, reflecting the 
fact that their total degree is nonpositive (an operator of this type is some- 
times referred to in the literature as exactly solvable since it preserves an 
infinite flag of polynomial subspaces). Note that the terms with coefficients 
ai7, 028 and 033 have positive y-degree, so they cannot be obtained from a 
polynomial in the generators of 0p^'^. This is one of the first examples of a 
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partially solvable differential operator preserving a finite-dimensional space 
of polynomials which is not Lie algebraic. (A number of partially solvable 
operators which preserve a finite-dimensional module of monomials but are 



not Lie algebraic are also discussed in [16|. The exponents of the monomials 



spanning these modules do not form a continuous chain {0, 1, . . . , n}). 



4 Matrix k-th order differential operators in 
with invariant modules of polynomials 

We shall now consider the more general case of A;-th order 2x2 matrix differ- 
ential operators admitting finite-dimensional invariant subspaces. Our goal 
is to prove structure theorems for these operators, which are the counter- 



parts of Theorems p.4| , |3.6| and p.lO| , in the case in which the invariant rank 2 
module C C°°(t/, C © C) is either a direct sum of rectangular, triangular 
or staircase modules. 

Consider the associative algebra 2) = of all 2 x 2 matrix differential 

operators with smooth coefficient functions. We introduce in 2) the usual 
Z2-grading; an operator 



Til Ti2 
T21 T22 



(4.1) 



is said to be even if T12 = T21 = 0, and odd if Tn = T22 = 0. The associative 
algebra D becomes then a Lie superalgebra with generalized Lie product 
given by 

[Ti, T2I = T1T2 - (_l)degTidegT2y^y^ _ (4 2) 

(We shall use the symbol 2) to denote both the associative algebra and 
the Lie superalgebra). We will construct graded subalgebras of D of matrix 
differential operators preserving the above mentioned vector- valued modules, 
in analogy with the scalar case. 



4.1 The rectangular module 7lni,mi ® ^^a ma 

This case has been recently studied in p|. We shall be very concise, for the 
results are a natural generalization of the single- variable case considered in |1| 
and [^. We shall assume that n = n2 > rii and m = mi > m2 (the other 
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choice 7712 > ""^1 leads to analogous results). Let A = n — rii and F = m — m2. 
Let us denote the direct sum 7ln-A,m © T^n,m-r by Consider the 

graded subalgebra s^;^ of 2) generated by the 8 + 2(A + l)(r + 1) matrix 
differential operators given by 



= diag( J^_^, J^) , = diag(ir^, K:,_r) , e = ±, , 

J = - diag(n + A,n) , K = - diag(m, m + F) , 

Qap = (la{x)qf3{y,m,T)a- , Ql^ = qa,{x,n, A) qp{y) a+ , (4.3) 

where a = 0, . . . , A , /3 = 0, . . . , F , 0"+ = (a^)* = f q g j > S'lid 



ga(a;) = , ga(x, n, A) = J]^ - n + A - /c) 9" , 

fc=i 

with the convention that a product with its lower limit greater than the 
upper one is defined to be 1. Let s^'^ be the abstract Lie superalgebra 
corresponding to the realization s^m (it '^^Y be shown that the supercom- 
mutators ( [4.2|) of the operators ( [4.3|) do not depend explicitly on n and m, 
|]1[). Just as in the single- variable case, 0, the Lie superalgebra s^'^ will 
be infinite-dimensional unless the gaps A and F are suitably restricted. The 
corresponding finite-dimensional Lie superalgebras are easily identified. 

Lemma 4.1 The Lie superalgebra s^'^ is finite- dimensional if and only if 
A + F < 1. We have: 

gho ^ ^0,1 ^ ^pj(2, 1) © , s°'° ~ © sli 

i=l 

Remarks. We have not included in s^'^ the central term given by K. We 
regard s°'° as a Lie algebra with J and K replaced by J = cxs. Finally, recall 
the isomorphism sp[(2, 1) ^ osp(2,2), [p^ . 



Let iH^ ^™ be the associative subalgebra of 2) spanned by the operators 

which preserve 7?.^^^. It is easy to see that s^'^ C ^^7m-r^ ^^^^ 
the action is irreducible. Thus, the Burnside Theorem guarantees that any 
matrix differential operator preserving this module may be expressed as a 
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polynomial in the generators of 5^'^ plus an operator annihilating the mod- 
ule. The kernel of the homomorphism '■ ^n~rn-r ~^ ^^^i'^rTm-?) 
simply the set of operators ( [4.1| ) whose component T^- belongs to keTHnj,mj 

as given in Lemma We choose a distinguished complement iH„ 
kerTT^;^ just as in the scalar case. The following lemma is a direct general- 
ization of the corresponding single- variable result, |^: 

Lemma 4.2 Let p^m '■ •^(s^''") — > 2) 6e the homomorphism determined by 
the representation p^'^ : s^'^ — >• s^^- The monomials 

{xis^Yis^nT^yiTy , Qt^ois^'YiT^yiTy , 

Ql.is^YisyiTy , Q'^,{s^y{Ty}, (4.4) 

where X = Qqq, 1, J, e = —e* = ±, a = 1, . . . , A, and P = 1, . . . ,T, form 
a basis o/imp^'^. 

Note that J may be expressed in terms of E, J, K (or for 3°'''). The 
number of monomials in (O) acting non-trivially in T^n'm-r precisely 

Theorem 4.3 Let T^^^ be a k-th order differential operator in ^nm-v- 
Then T^''^ may be expressed as a linear combination of the monomials in ( [4.4|) 
of differential order at most k. 

4.2 The triangular module 7^^ © 7^^ 

We shall assume that n = n2 > fii, and let A = n — rii. Let us denote the 
module T^-a © by Let T^~^ be the associative subalgebra of 2) 

of differential operators preserving r^^^. Consider the 9 + (A -M)(A -f 2) 
matrix differential operators given by 

T = diag( j;) , J = i diag(n + 2A,n), z = 1, . . . , 8 , 
Qap = 1af3(^~, Qtt3 = qapin, A) a+ , < « + /3 < A , (4.5) 
with given in (|2.3|), and 



A-a-/3 

g„/j = a;V, ga/3(n,A)= J] (x^, + - n + A - A;) 9^ . 

fe=i 
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We observe that the even generators in ( [4.5| ) span a Lie algebra which is 
isomorphic to qI^. It may be easily verified that both the generators Q~p 
and span irreducible modules under the adjoint action of the even sub- 
algebra, which are isomorphic to the standard cyclic s la-module of highest 
weight (A,0), with J acting as a scalar. Moreover, the anticommutator of 
and Q^^ can be expressed as a polynomial of degree A in the even gen- 
erators and the identity, depending explicitly on n only through a multiple 
of the identity. Let us denote by the Lie superalgebra generated by the 
differential operators ( [4.51 ). 

Lemma 4.4 The Lie superalgebra is finite- dimensional if and only if 
A < 1. We have: 

4~sp[(3,l), S°~S[3©S[2. 

Remark. We regard s° as a Lie algebra with J replaced by J = (T3. 

In general, the Lie superalgebra of differential operators preserving a pair 
of A^- dimensional "triangular" modules with a gap A = 1 is isomorphic to 
spi{N + l, 1). It is easy to see that C X^"^, and that the action on T"~^ 
is irreducible. The kernel of the homomorphism ir^ : — > End(7^"~^) 
is just the set of the operators ( |4.1| ) with entries Tjj in kervTn. as given in 
Lemma p.2| . As usual, we choose as distinguished complement to ker vr^ the 

vector space X„ spanned by the differential operators in with com- 

ponents of the form Tij = Es+t<nj fst{x, y) <9^9*. 

We first note that the definition of bidegree can be naturally extended to 
matrix- valued differential operators acting on C[x, ©C[x, y]. The bidegrees 
of the generators ( |4.5| ) of are well-defined. The bidegree of each is 
given by the bidegree of the corresponding J^, and deg J = (0,0). For the 
generators of the odd subspace, we have degQ^^ = ^{a,(3). In order to 
define the length of a monomial in the generators (|4.5| ) we need to be a little 
more careful. Let = (T^)^^ ■ ■ ■ (T^)'^**. Consider the monomials of the 
form Tf = X, T^, where e = 1, 2, ±, X, = Pi, = P2, X± = and Pi, 
P2 are the canonical projectors, which can of course be expressed in terms of 
II and J (or J, if A = 0). We define the length of a monomial as 

\T^\ = min \L\ , 

degT/'=degT/f 

where the minimum is taken over the set of monomials with \L\ < \K\. 
Again, a monomial is of maximal length if \T^\ = \K\. Thus, Pi, P2 and 
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the Qaf3^s are all of maximal length (equal to zero). For e = ±, the bidegree 
does not determine uniquely the monomials of maximal length. The structure 

• • • • 

• o o o • ® Length I = 

• o Length I = 1 

• • Length / = 2 
• o o • 



® 



® 

(0,0) 

® 



Figure 3: Lattice of bidegrees of monomials T_p of length / < 2 for A = 1. 

- n—A 

theorem for may now be stated as follows: 

Theorem 4.5 Let T^''^ be a k-th order differential operator in X„ . Then 
T^^'^ may he expressed as a linear combination of the monomials 

{T^ {T'YiTy \e=l,2,±, 0<i+j + \L\<t<k- , (4.6) 

where are monomials of maximal length and different bidegrees for each 
(fixed) value of e. 

Remark. If A > k, no monomial {T'^y{T^y may appear in the expansion 
of TW. 

Proof. The monomials ( [4.6|) are linearly independent. Let denote the 
number of monomials of order at most k and type e in ( ^4.6| ). We already 
know that N^''^ = \{k + ifik + 2)^, cf. (p]). Let r(no) = dimT^^, = 
i(no + l){nQ + 2). An easy graphical argument, cf. Fig. 3, shows that 

N^=^-^^^^^r{k-6,^A)+ ^ 3{2l + A)r{k-6,^A-l) 

1=1 

= ^{k + l){k + 2){k tA + l){k T A + 2) . (4.7) 

— n—A 

Therefore, N^_^ + + + N~_^ = dim X„ . The theorem follows from 
the fact that the highest order derivatives of the monomials (|4.6| ) are all in- 
dependent. Q.E.D. 
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The dimension of the subspace ' of of differential operators 

of order at most k is given by 

dimtr^'^''^ = ^{k +l){k + 2){2e + Qk + 4 + A^) . 

It is worth mentioning that in the one- dimensional rank two case this number 
is independent of A, 



4.3 The staircase module S'L „ (B SI. ^ 

For the sake of simplicity, we shall restrict ourselves to the case p = pi + 1 = 
P2, 1 = 11 = 12- Let us denote the module ^ © S^^^ by S'p^q- Consider 
the finite-dimensional Lie superalgebra spanned by the 10 + 2r matrix 
differential operators given by: 

T = diag( j;) , J = ^ diag(p + 1, p) , ^ = 1, . . . , 5 + r , 

Qq = a' , Qi = xa' , Qo = i^d^ + rydy - p) , Qt = dx cr^ , 

Qt+o. = ^''dy^^ . a = 0,...,r-l, (4.8) 



with Jp given in ( |2.4| ). The abstract structure of the Lie superalgebras 
is the semidirect product pl(2, 1) k C'"^^''", where C^"*"^'^ corresponds to the 
Abehan Lie superalgebra spanned by T^+* and Qt+a^ z = 0, ... ,r, a = 
0, . . . , r — 1. It may be readily verified that is contained in the associative 
subalgebra ^ of 2) of differential operators preserving 5^ ^. The kernel 
of the homomorphism VTp'^J : — > End(«Sp^^) is spanned by the differential 
operators (|4.1| ) with components Tij in ker TTp, as given in Lemma p.2| . We 



~ r 1 

choose as a distinguished complement to kervr^'g the subspace ©^'^ C &p^q 
of differential operators with components of the form: 



0<s+rt<pj 
0<t<q 

We note that just as in the scalar case, the module «Sp ^ is reducible (but not 

- r,l 

completely reducible) under the action of Sp. In fact, no operator in &p^q 
with positive ^/-degree may be obtained from a polynomial in the generators 

* r, 1^ " r,l 

of Sp. Let ©p g denote the subspace of ©^ ^ of differential operators with 
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nonpositive y-degree. The structure theorem for &p^q may be proved by 
essentially the same argument used in Theorem p.lO| . 

Theorem 4.6 Let T*-'^^ be a k-th order differential operator in &pq ■ If 
k<q, then T^^^ may he represented as a {k + l)-th degree polynomial in the 
generators of Sp. If k > q, then T*^*^) may be expressed as the projection of 

~ r 1 

such a polynomial in &pg along kervr^'g. 



Proof. The key of the proof is to obtain a suitable set of monomials in 

~ r 1 

the generators of Sp such that their projections in ©p'^ along ker tt^'^ are all 
independent. We consider monomials of the form 

rpL {T^yijy^ (4 9) 

where e = 1, 2±, and are monomials of maximal length and different bide- 
grees for any (fixed) value of e. For e = 1,2, Lemma |3.9| implies that we can 
obtain dim ©pf^g monomials of the form ( |4.9| ) with independent projections in 

&p^q along ker Tr^'^^. Let &p^_q denote the subset of &p^q of differential oper- 
ators mapping Pi(«Sp'^) into P2(«5p g). The projections in &p q along keryr^'^^ 
of the monomials 



(4.10) 

with e = 0, 1, < J < g, < s < p - 1 - jr, < t < J - 1, 1 < i < r - 5to, 
0<n,0<m<q — j, and n + rm < p — s — jr (where s = s for the first 
group of monomials, and s = for the second group), form a basis of ©pfq~- 
The number of monomials in ( [4.10| ) is given by 

g p-i-j> _ ^ 
Kg =^ Yl s)(T{r,p-l - s - jr, min (g, [ ] ) - j) , 

j=0 s=0 ^ 

where f{j, s) = 2 if s > 1 and /(j, 0) = jr + 2, and a{ro,po, go) = dimS^lq^. 
It may be verified that 

iV;- = dim = ^^1 + ^^^^ + ^^ (^p(p + 1) _ qr{p + ^) + ^(3g + 1)) . 
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We also leave as an exercise to the reader to obtain 

at;; = dim e;f; = + ^^^^^ + {p{p + l) - gr(p + ^) + ^(3g + 1)) 

— r 1 

monomials of the form (T'^)*(T'^)-' with independent projections in &pq 
along kervTp'^. The theorem follows from the fact that the highest order 
derivatives of the monomials ( [4 .91 ) are all different. Q.E.D. 
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